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We propose the design of a space-variant Wien filter for electron beams that induces a spin 
half-turn and converts the corresponding spin angular momentum variation into orbital angular 
momentum (OAM) of the beam itself, by exploiting a geometrical phase arising in the spin ma- 
nipulation. When applied to a spatially-coherent input spin-polarized electron beam, such device 
can generate an electron vortex beam, carrying OAM. When applied to an unpolarized input beam, 
the proposed device in combination with a suitable diffraction element can act as a very effective 
spin-polarization filter. The same approach can be also applied to neutron or atom beams. 



Phase vortices in electronic quantum states have been 
widely investigated in condensed matter, for example in 
connection with superconductivity, the Hall effect, etc. 
Only very recently, however, free-space electron beams 
exhibiting controlled phase vortices have been experi- 
mentally generated in transmission electron microscope 
(TEM) systems, using either a spiral phase plate ob- 
tained from a stack of graphite thin films [l[ , or a "pitch- 
fork" hologram manufactured by ion beam lithography 
[H, 01] . In a cylindrical coordinate system r, 4>, z with 
the z axis along the beam axis, a vortex electron beam 
is described by a wavefunction having the general form 
4>(r, <j), z,t) = u{r, z,t) exp(i£<j)), where I is a (nonzero) 
integer and u vanishes at r = 0. As in the case of atomic 
orbitals, I is the eigenvalue of the ^-component orbital 
angular momentum (OAM) operator L z — — id^, (in units 
of the reduced Planck constant Ti) , and therefore an elec- 
tron beam of this form carries Ifi of OAM per electron Q ■ 
The recent experiments on electron vortex beams were in- 
spired by the singular optics field, in which similar phase 
or holographic tools have been used in the last twenty 
years (see, e.g., jl| and references therein). In optics, a 
recently introduced alternative approach to the genera- 
tion of vortex beams is based on the "conversion" of the 
angular momentum variation occurring in a spin-flip pro- 
cess into orbital angular momentum of the light beam, 
when the latter is propagating through a suitable spa- 
tially variant birefringent plate [1,0|- In this paper, we 
propose that a beam of electrons traveling in free space 
undergoes a similar "spin-to-orbital angular momentum 
conversion" (STOC) process in the presence of a suitable 
space- variant magnetic field. The same approach may 
work also for neutrons, or any other particle endowed 
with a spin magnetic moment (e.g., atoms, ions). Of 
course, in the case of electrons, as for other charged par- 
ticles, the magnetic field, besides acting on the spin, will 
also induce forces that must be compensated in order to 
avoid strong beam distortions or deflections. Such com- 
pensation may be obtained by a suitable electric field, 



and this leads us to conceiving the proposed apparatus 
essentially as a space-variant Wien filter. Such appara- 
tus can be exploited for generating vortex electron beams 
when a spin-polarized beam is used as input. Conversely, 
if a pure vortex beam is used in input, by means e.g. of a 
holographic method, one can use the STOC process for 
filtering a single spin-polarized component of the input 
beam, as we will show further below. 

Let us consider an electron beam propagating in vac- 
uum along the z-axis and crossing a region of space lying 
between z = and z = L in which it is subject to electric 
and magnetic fields E = — V<I> and B — V x A, where 
4> and A are the scalar and vector potentials, respec- 
tively. In the non-relativistic approximation and neglect- 
ing all Coulomb self-interaction effects (small charge den- 
sity limit), the electron beam quantum propagation and 
spin evolution are generally described by Pauli's equation 
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where ip is the spinorial two-component wave-function 
of the electron beam, e = — |e| and m are the electron 
charge and mass, dt is the derivative with respect to the 
time variable t, ft, = — ^gustr is the electron magnetic 
moment, with hb = h\e\/2m the Bohr's magneton, g ~ 
2 the electron g-factor, and a = {(T x ,d y ,d z ) the Pauli 
matrix vector. 

As a first step, we consider the simpler case in which 
the electric and magnetic fields are taken to be uniform, 
lying in the transverse plane xy, and arranged as in stan- 
dard Wien filters 0,0], be. perpendicular to each other 
and balanced so as to cancel the average Lorentz force, 
i.e. Eq — Bop c /m where Eg and Bq are the electric and 
magnetic field moduli, and p c the average beam momen- 
tum. The magnetic field B is also taken to form an arbi- 
trary angle a with the axis x within the xy plane. For this 
case, we solved the full Pauli's equation in the paraxial 
slow-varying-envelope approximation for an input beam 
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having a gaussian profile and an arbitrary uniform input 
spin state |V>)i« = a i|t)+ a 2|l), where ||) and \l) denote 
a state for which the spin is parallel or antiparallel to the 
z axis, respectively. The complete expression of the re- 
sulting spinorial wave- function is given in the supplemen- 
tal material (SM) [l(|, while here we summarize the main 
findings. The beam propagation behavior corresponds to 
the well known astigmatic lensing in the plane perpen- 
dicular to the magnetic field. More precisely, the beam 
undergoes periodic width oscillations, with a spatial pe- 
riod A2 = ttR c , where R c = p c /(\e\Bo) is the cyclotron 
radius. This lensing phenomenon is also predicted by a 
classical ray theory, when properly taking into account 
the effect of the input fringe fields [lOj . The output spin 
state is instead given by the following general expression 
(Eq. 4 in SM 0) 

\ip)out = «i [cos(«5/2)|t> +sin(<S/2)ie iQ ||>] 

+a 2 [cos(V2) I i) + sm(5/2)ie- ta | f)] , (2) 

where S — AttL/Ai and Ai = 4irR c /g ~ 2A2. This spino- 
rial evolution corresponds to the classical Larmor preces- 
sion of the spin with spatial period Ai/2, 5 being the 
total precession angle. However, in addition to the spin 
precession, Eq. ([2]) predicts the occurrence of wavefunc- 
tion phase shifts. In particular, for a | f) or | |) input 
state and a total spin precession of exactly half a turn, 
i.e. 8 = 7T or L = Ai/4, the wavefunction acquires a 
phase shift given by ±a + 7r/2, where a is the magnetic 
field orientation angle mentioned above and the ± sign 
is fixed by the input spin orientation (+ for | f) and — 
for These phase shifts can be interpreted as a spe- 
cial case of geometric Berry phases arising from the spin 
manipulation (llj . 

Let us now move on to the case of a spatially vari- 
ant magnetic field. We consider multipolar transverse 
field geometries with cylindrical symmetry, described by 
the following expression for the magnetic field (with 
the vector given in cartesian components): B(r,(f),z) — 
Bo(r) (cos a(4>), sin a(0),O), where the angle a is now the 
following function of the azimuthal angle: 



a(r,<f>,z) =q(j> + l3 
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where q is an integer and /3 a constant. Clearly, such 
a field pattern must have a singularity of topological 
charge q at r — 0. In particular, by imposing the vanish- 
ing of the field divergence, we find that the radial factor 
Bo(r) ~ r~ q , i.e., the field vanishes on the axis for q < 0, 
while it diverges for q > 0. In the latter case, there must 
be a field source on the axis. We call "g-filter" a balanced 
Wicn filter whose magnetic field distribution in the beam 
transverse plane obeys Eq. (|14l) . The electric field will be 
taken to have an identical pattern, except for a local 7r/2 
rotation, so as to balance the Lorentz force. Some exam- 
ples of such g-filter field distributions are shown in Fig. 
[TJ We are particularly interested in the negative q ge- 




FIG. 1. (color online) Electric (upper panels) and magnetic 
(lower panels) field q-filter geometries for different topological 
charges: (a) q = —2, (b) q = —1, (c) 9=1, and (d) q = +2; 
in all cases /3 — tt/2. 



ometries, which do not require to have a field source at 
r = 0. For example, the q = — 1 case corresponds to the 
standard quadrupole geometry of electron optics, while 
q = — 2 corresponds to the hexapole one. Wien filters 
with such geometries have been already developed in the 
jast for the purpose of correcting chromatic aberrations 
T2I 13j . Moreover, inhomogeneous Wien filters including 
several multipolar terms have been also considered for the 
purpose of spin manipulation, with the added advantage 
of obtaining a stigmatic lensing behavior 14{ . A possible 
design of the q = — 1 filter with quadrupolar geometry is 
shown in Fig. [2] In such non-uniform field geometry we 
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FIG. 2. (color online) Electrodes and magnetic poles geom- 
etry of a g-filter with q = — 1 (quadrupole), seen in cross- 
section (a) and in three-dimensional rendering (b). The filter 
length is set to 50 cm. In panel (a), the calculated vector 
potential A z (in false colors), which is also roughly propor- 
tional to the electric potential, and the projection on the xy 
plane of the simulated electron ray trajectories for 100 keV 
energy are also shown (see SM for details about the simula- 
tions (ioj]h for a ring-shaped input beam with radius r = 100 
fim (ray color grows darker for increasing z). In panel (c) 
a zoomed-in view of the central region. The magnetic field 
at r needed to obtain the tuning condition 8 = -n is 3.5 mT, 
with a corresponding electric field of 575 kV/m. These are 
obtained with an electrode potential difference of « 9 kV and 
a magnetization of 135 A/mm. The fields need to be set to 
the design- values with a precision of 1 part in 10 4 . 
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cannot solve analytically the full Pauli's equation. How- 
ever, the beam propagation is already well described by 
classical dynamics and can be derived either analytically, 

or by numerical ray 
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using a power-expansion m r 
tracing. In the former case, we find that to first order the 
q- filter for q ^ is already stigmatic, i.e., it preserves the 
beam circular symmetry. Only second-order corrections 
introduce aberration effects [141] . This behavior is con- 
firmed by our numerical ray-tracing simulations (see SM 
for details [lO]), which show relatively weak higher-order 
aberrations (see Fig. [5]). It should be noted that these 
simulations have been performed for realistic values of the 
electric and magnetic fields, as required to obtain a spin 
precession of half a turn across a propagation distance of 
50 cm at a beam radius of 100 /im. These calculations 
are expected to reproduce very well the electron density 
behavior (and spin precession) as would be obtained from 
Pauli's equation. However, Pauli's equation predicts an 
additional purely quantum phenomenon, namely the geo- 
metric phase already discussed above. In a semi-classical 
approximation, the geometric phase will still be given by 
±a + 7r/2, where a is however now position-dependent 
and given by Eq. (fT4|) . More specifically, neglecting the 
aberrations, each possible electron trajectory within a 
beam is straight and parallel to the z axis. Therefore, 
the electrons travelling in a given trajectory will expe- 
rience a constant magnetic field of modulus B^{r) and 
orientation a(cf>). The set of all electrons travelling at 
a given radius r will then undergo a uniform spin pre- 
cession by angle 5(r) and, if 5(r) = it (i.e., for a spin 
half-turn rotation), they will also acquire a space-variant 
geometric phase given by ±a(</>) + tt/2 = ±q<p± /3 + 7r/2, 
with the ± sign determined by the input spin orienta- 
tion. In other words, the outgoing wavefunction acquires 
a phase factor exp(i£</>), with I — ±q, corresponding to a 
vortex beam with OAM ±qh. In a quantum mechanical 
notation, spin-polarized input electrons with given initial 
OAM £ passing through a g-filtcr undergo the following 
transformations: 

| tJ) -> cos(<5/2)| t, £) + ie ifs sin(5/2) \±,£ + q) 
\i,£) -> cos(5/2)\l,£)+ie- ll3 sm(S/2)\t,£-q) (4) 

where the ket indices now specify both the spin state 
(arrows) and the OAM eigenvalue. 

Equations (@| show that, in passing through the q- 
filter, a fraction / = sin 2 (<5/2) of the electrons in the 
beam will flip their spin and acquire an OAM ±qh, while 
the remaining fraction 1— / = cos 2 (£/2) will pass through 
the filter with no change. When L = Ai/4, then S = n 
and all electrons are spin-flipped and acquire the cor- 
responding OAM. In the specific case q = 1, the spin 
angular momentum variation for the electrons undergo- 
ing the spin inversion is exactly balanced by the OAM 
variation, so that the total electron angular momentum 
remains unchanged in crossing the filter. This is the pure 
"spin-to-orbital conversion" STOC process mentioned in 



the introduction, and it occurs for q = 1 because this 
geometry is rotationally invariant and therefore no angu- 
lar momentum can be exchanged with the field sources 
in the filter. In the g ^ 1 case, the input spin still con- 
trols the sign of the OAM variation, but the total beam 
angular momentum is not conserved and some angular 
momentum is exchanged with the field sources. We note 
that this OAM variation can be also explained as the 
effect of the spin-related magnetic-dipole force acting on 
the electrons within the magnetic field gradients, as more 
fully discussed in the SM [10|. 

The "tuning" condition L = Ai/4or<5 = 7r can be 
achieved in principle for a given radius r by adjusting the 
strength of the magnetic and electric fields or the device 
length L. Since the precession angle S is r-dependent, 
however, this tuning condition can be applied to the en- 
tire beam only if it is shaped as a ring, i.e. with all elec- 
tron density peaked at a given radius r. Vortex beams 
with OAM I 7^ typically have a doughnut shape, so 
they approximate a ring fairly well. On the other hand a 
gaussian input beam (with i = 0) cannot be fully trans- 
formed, as S = at r = 0, where the beam has the 
maximum density. In such cases, only a fraction / of the 
electrons would be converted. 

So far we have assumed a spin-polarized input beam. 
However, high brightness (i.e., spatially coherent) spin- 
polarized electron beams, suitable for high-resolution 
TEM applications, are not so easily available. State-of- 
the-art spin-polarized sources may achieve a brightness 
of 10 7 A cm -2 sr _1 and a polarization purity of up to 
90% [IH (and the source decays with time due to laser- 
induced damage). It is interesting then to analyze the 
effect of the g-filter on a initially unpolarized electron 
beam, having arbitrary initial OAM £. Such input can 
be simply viewed as a statistical mixture in which 50% 
of the electrons are in the state | f, ^) an d 50% in the 
state | After passing through a tuned g-filter, the 

beam becomes a 50-50 mixture of states | I, £ + q) and 
| f, £ — q), for which spin and OAM are correlated (if the 
(/-filter is not tuned, the fraction of converted electrons 
decreases to f/2 in each spinorbit state, and there will be 
a residual 1 — / fraction of electrons in states | f, £) and 
|4-,^))- As we discuss now, this spin-OAM correlation can 
be exploited for making an effective electron beam spin- 
polarization filter. Such filter requires four basic elements 
in sequence (as shown in Fig. 2 of SM [l3|): (i) an OAM 
manipulation device, such as a fork hologram 0,(3], to 
set £ =^ 0; (ii) a g-filter with q = £, generating a mixture 
of electrons in states \i,2£) and |f, 0); (hi) a free propa- 
gation (or imaging) stage that allows these two states to 
develop different radial profiles by diffraction, because of 
their different OAM values; (iv) a circular aperture for 
finally separating the two states. In particular, in stage 
(iii) state 1 4-, 2£) will acquire a radial doughnut distribu- 
tion, as in Laguerre-Gaussian modes with OAM 2£, which 
vanishes close to the beam axis as r 2£ , while state 1 "f, 0) 



4 



will become approximately gaussian, with maximum in- 
tensity at the beam axis, as shown in Figs. [3^i,b. There- 
fore, a suitable iris (Fig. [3};) will select preferentially the 
electrons in the fully polarized state | y, 0) . An opti- 
cal OAM sorter exploi ting a similar approach has been 
demonstrated recently [16| . A specific calculation for the 
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FIG. 3. (color online) Electron beam profiles in the far field 
of the g-filter for the I = component (panel a) and the 
I = 2 component (panel b), and a possible discriminating iris 
radius (panel c) r = w to be used to separate them in order to 
make a spin-polarization filter. Panel (d) shows the intensity 
profiles of the same components (t = - blue dot-dashed 
line, 1 = 2- green dashed line) and of the possible residual 
I = 1 component for an untuned g-filter (red solid line), w 
is the gaussian beam waist radius in the far-field plane. A 
realistic value for the iris radius is of the order of several tens 
of microns (obtained by setting the aperture some distance 
after the focal plane of the second condensor) 

case \q\ = \i\ = 1 and an iris radius equal to the beam 
waist w in the "far field" yields a transmission efficiency 
of our device of 55.5% (not including the losses arising in 
the OAM-manipulation device) and a polarization degree 
(Jf — Ii)/(If + Ii), where are the two spin-polarized 
currents, of ~ 97.5%. Higher degrees of polarization can 
be obtained at the expense of the efficiency by reducing 
the iris diameter or by employing higher q values (or vice 
versa). It is worth noting that this apparatus works also 
with a partially tuned g-filter, as in this case the unmod- 
ified electron beam component is left in the initial OAM 
state £ = q and therefore is also cut-away by the iris. An 
untuned g-filter will however have an efficiency reduced 
by the factor / = sin 2 (<5/2). The aberrations introduced 
by the g-filter, even if left uncorrected, might also affect 
its efficiency but not its main working principle, as this 
is based on the vortex effect, which is protected by topo- 
logical stability. Finally, the possible spin depolarization 



effect of fringe fields can be neglected if the length-to-gap 
ratio of the filter is large enough [IT 
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We note that the spin-filter application discussed above 
is a new counterexample of the old statement by Bohr 
that free electrons cannot be spin-polarized by exploi ting 
magnetic fields, due to quantum uncertainty effects [171 - 
fig( | . The reason why we can overcome Bohr's arguments 
is essentially that we do not use the magnetic forces di- 
rectly to obtain the separation, but take advantage of 
quantum diffraction itself, as also proposed recently in 
Ref. [13 (see SM for a fuller discussion 0). 

In conclusion, we believe that the g-filter device de- 
scribed in this paper can be manufactured relatively sim- 
ply, for applications in standard electron beam sources 
such as those used in TEMs or other kinds of electron mi- 
croscopes. In combination with current field-effect unpo- 
larized electron sources, such filter might provide a spin- 
polarized source with a brightness ~ 10 9 A cm~ 2 sr , 
about two orders of magnitude higher than the current 
state of the art. This result, if it will be proved practical 
enough, may open the way to a spin-sensitive atomic- 
scale TEM, e.g. suitable for investigating complex mag- 
netic order in matter or for spintronic applications. 

We acknowledge the support of the FET-Open Pro- 
gram within the 7 th Framework Programme of the Euro- 
pean Commission under Grant No. 255914, Phorbitech. 



APPENDIX: SUPPLEMENTARY MATERIAL 

Gaussian electron beam in homogeneous Wien filter 

Let us consider an electron beam propagating in vac- 
uum along the z-axis, subject to transverse and or- 
thogonal electric and magnetic fields E and B lying 
in the xy plane. In this Section, we assume the two 
fields to be spatially uniform, as in standard Wien fil- 
ters, but we will release this assumption later on. If 
the magnetic field direction makes an angle a with 
the x-axis, we may write E = i?o(sin a, — cos a, 0) and 
B = £?o(cos a, sin a, 0). As scalar and vector poten- 
tials, we may then take $ = —Eq(x sina — ycosa) and 
A = Bq (0,0,ycosa — a; sin a), respectively. In the non- 
relativistic approximation and neglecting all Coulomb 
self-interaction effects (small charge density limit), the 
electron beam propagation is described by Pauli's equa- 
tion 



ihdttp = 



-L(- t hV-eA) 2 
Zm 



e<I> - B ■ p, 



4' (5) 



where ip is the spinorial two-component wave-function 
of the electron beam, e = — j e j and m are the electron 
charge and mass, dt is the derivative with respect to the 
time variable t, ft — —\giiBG is the electron magnetic 
moment, with \ib = h\e\/2m the Bohr's magneton, g ~ 
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2 the electron g-factor, and <r = (a x , fry, fr z ) the Pauli 
matrix vector. 

We seek a monochromatic paraxial-wave solution 
with average linear momentum p c and average energy 
E c = p 2 ./2m in the form ip(x,y, z,t) — exp[i(p c z — 
E c t)/h]u(x,y, z), where u(x, y,z) is taken to be a slow- 
varying envelope spinor field. Inserting this ansatz in 
Eq. ([5]) and neglecting small terms in the z derivatives of 
u, we obtain the paraxial Pauli equation 



In addition to the spin dynamics, from Eq. (fT2)) . we 
find that the wavefront complex curvature radius q 2 (z) 
in the direction perpendicular to the magnetic field B 
changes periodically with spatial period A2 . These oscil- 
lations correspond to the well known astigmatic lensing 
effect of the Wicn filter. 



Space-variant Wien filter and ray-tracing simulations 



(6) 



where = d 2 



dy is the Laplacian in the beam trans- 



verse plane and k c = p c /h is the De Broglie wavevec- 
tor along the beam axis. In deriving Eq. ([6]) we set 
Eq = pcBo/m, so as to cancel out the average Lorentz 
force on the beam. 

Equation ([SJ) determines the evolution of the spinor 
field u along the axis z. This equation can be solved 
analytically for the case of a gaussian input beam, as 
given by u(r, <fi, 0) = aexp[—r 2 /wq] for z = 0, where 
Wo is the beam waist and a = (01,02) the input spinor, 
corresponding to an arbitrary spin superposition \ip) = 
Oi| t) + 02I D- A straightforward calculation shows that 
the solution is given by 



u(r, 4>, z) — G(r, </>, z)M(z)a 
where M(z) is the matrix 
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where Ai = (27rh 2 k c )/ (mgfisBo) and the scalar gaussian 
factor G{r, <p, z) is given by 



G{r,(j),z) = g (z)e 



1 2,i(z) + 2„ 2 (z) 
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fl Ai/4~ Ai/2. 



The matrix M(z) given in Eq. ([8]) corresponds to Eq. 
2 in the main paper and describes the magnctic-ficld- 
induced Larmor spin precession during propagation and 
the associated geometric phases discussed in the main 
paper. The precession length Ai corresponds to two full 
spin rotations. 



In this Section, we consider "g-filter" field geometries 
with cylindrical symmetry, described by the following ex- 
pression for the magnetic field (with the vector given in 
cartesian components) : 

B{r,(j),z) = B Q {r)(co$a{<f>),sm.a(4>),Q) (13) 

where the angle a is now the following function of the 
azimuthal angle: 



a(r,<f>,z) =q4> + fi 



(14) 



where q is an integer and /3 a constant. As discussed in 
the main paper, such a field pattern must have a singu- 
larity of topological charge q at the beam axis r = 0. In 
particular, by imposing the vanishing of the field diver- 
gence, we find that the radial factor Bo(r) ~ r~ q , i.e., the 
field vanishes on the axis for negative values of q, while 
it diverges for positive values of q. In the latter case, 
there must be a field source on the axis. The electric 
field will be taken to have an identical pattern, except 
for a local tt/2 rotation, so as to balance everywhere the 
Lorentz force. We are particularly interested in the neg- 
ative q geometries, which do not require to have a field 
source at the beam axis. For example, the q = — 1 case 
corresponds to the standard quadrupole geometry, while 
q = — 2 corresponds to an hexapole one. 

The electron propagation problem in such non- uniform 
field can only be solved approximately or numerically. A 
second-order geometrical optics solution is reported in 
Ref. [? ] , which includes also a detailed analysis of spin 
precession but without considering the geometric phase 
effect (as well as the magnetic dipole forces). This anal- 
ysis shows that, to first order, a balanced Wien filter 
having quadrupole or hexapole geometry but vanishing 
dipole term (as in our g-filter) is already stigmatic, i.e. 
it preserves the cylindrical symmetry without need for 
further corrections. The only beam lensing or distortion 
effects enter as higher-order aberrations (i.e., the G term 
in Eq. (31) of Ref. [? ]). 

To further analyze the behavior of our g-filter, we have 
performed ray-tracing simulations of the electron propa- 
gation for a quadrupole geometry (q = —1). The mag- 
netic and electric fields have been calculated using COM- 
SOL finite-elements simulation (www.comsol.com). The 
ray-tracing routines therein contained have been used to 
assess the shape of the beam at the exit of the filter start- 
ing from an input beam having 100 keV of enery that is 



6 



shaped as a ring, with a radius r of 100 /im. The magnetic 
field at r needed to obtain the tuning condition S = tt is 
3.5 mT, with a corresponding electric field of 575 kV/m. 
These are obtained in our geometry with an electrode po- 
tential difference of w 9 kV and a magnetization of 135 
A/mm. The fields need to be set to the design- values 
with a precision of 1 part in 10 4 . The electron velocity 
to be used as initial conditions in the simulations must 
be obtained from a model of the fringe fields. For the 
simplest description of the input fringe fields, i.e. the 
so-called sharp cut-off fringing field (SCOFF) model, the 
transverse velocity components remain constant and nil, 
i.e. v x — v y — 0, while the z-component v z (x,y) "just 
inside" the filter is given by the following expression: 

e$(x,y) _ eA z (x,y 



Vz(x,y) = v ~ 



mvo 



vo 



(15) 



where Vq is the velocity outside the filter. This condi- 
tion is derived from the conservation of energy and also 
ensures the conservation of the canonical momentum of 
the electrons entering the Wien filter and hence the con- 
servation of their wavefront orientation. The results of 
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FIG. 4. Projection on the xy plane of the simulated trajec- 
tories of the electrons starting from a circle of radius r = 100 
fim. The trajectories are shown with a color growing darker 
as z is increased. 

ray-tracing simulations are reported in Fig. 2 of the main 
paper and in Fig. [5] of the present manuscript, which 
show how a circle of electron positions representative 
of the input beam becomes deformed during the prop- 
agation into a quadrupole-lobed shape, due to second- 
or higher-order aberrations. For optimizing the perfor- 
mances of our proposed device, in particular in the spin- 
polarization filter application, these aberrations should 
probably be compensated by additional electron optics 
(e.g., octupole electrostatic lenses). Nevertheless, even 
for no aberration compensation, the vortex effect giving 
rise to the radial separation of the two spin components 
discussed in the spin-polarization application is expected 
to be preserved, owing to its topological stability. 



Magnetic dipole forces and angular momentum 

We now evaluate the effects of the force associated with 
the magnetic dipole of the electrons in the magnetic field 
gradient of the Wien filter. This is the force used in 
Stern-Gerlach experiments and the only spin-dependent 
force acting on the electrons, so it is important to assess 
its effects in detail. At any given point inside the filter, 
assuming that the magnetic dipole of the electron is /x = 
—gfisS, where S is the corresponding spin vector, the 
associated magnetic force is 



F = V(/x • B) 



(16) 



where the gradient must be taken while keeping fj, con- 
stant, i.e. acting only on the magnetic field coordinates. 
For our classical estimate, the magnetic dipole evolution 
can be described by the precession equation 



~dt 



9\e[ 
2m 



B x fi 



(17) 



Assuming a starting magnetic moment parallel to the 
z axis, the magnetic dipole will rotate around the local 
magnetic field direction. In the approximation in which 
aberrations are neglected, each electron travels along a 
parallel ray with constant transverse coordinates, and 
hence sees a uniform magnetic field, given by Eq. (|13p . 
Therefore, the magnetic dipole at any given point inside 
the filter will be given by the following expression: 

H(r, 4>,z) = ijlb[— sin<5(r, z) sin a((f>) , sin S(r, z) cosa(0), 
cosS(r,z)] (18) 

where S(r, z) is the total spin precession angle at z, for 
a trajectory at radius r. Taking the scalar product be- 
tween fi and B and then the gradient with respect to the 
coordinates of B only, we find that the magnetic dipole 
force is only directed along the azimuthal direction and 
given by 



2 ^sSo(r) . 
F = cpq — sin o[r, z) 



(19) 



where <j> denotes the unit vector along the azimuthal di- 
rection and we have used Eq. (|14|) for taking the <f> deriva- 
tive. 

Equation (fl~9| can be used for two purposes. First, we 
can verify that the ray deflection inside the Wien filter 
due to this force is negligible. Indeed, integrating the 
force in the time T needed for a full n rotation of the spin 
(or spin flip), we obtain the overall momentum change 



F+{t)dt 



2m 



g\e\B, 



Mr 



dS 



(20) 
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We can see from this equation that the momentum vari- 
ation is of the same order as the quantum uncertainty 
in momentum of the electron wavefunction, taking into 
account that the position uncertainty is of the same or- 
der of the beam size r (as it is also expected based on 
Bohr impossibility theorem, as discussed further below). 
This momentum change produces negligible effects in the 
electron propagation across the filter, although it instead 
becomes relevant in the far field. Indeed, it is precisely 
this force that produces the orbital angular momentum 
change at the basis of our spin selection. Indeed, we will 
now show that the torque associated with the magnetic 
dipole force is that responsible for the OAM change of 
the electrons. The z-component torque of this force is 
given by 



M z = (r x F) z = qn B B (r) sin 5 = q(fi x B) z 
dS z 



= -q- 



dt 



(21) 



where S z is the z component of the electron spin angular 
momentum S. Since this torque acts on the OAM of the 
electrons L z , we find the following angular momentum 
coupling law: 



dLz 
dt 



-Q- 



dS* 
dt 



(22) 



Once integrated for the entire flight time, corresponding 
to a full spin flip, this coupling law returns the same 
result reported in the main paper for the OAM variation 
induced by the spin flip. 

This correspondence proves that the effect, and the 
only effect, of the magnetic dipole force acting on the 
electrons is the change of OAM already considered by 
including the geometric phase shift in the wavefunction 
emerging from the filter. Therefore, the only significant 
effect of this force appears in the far field and is exactly 
the effect that we propose to utilize to separate electrons 
according to their spin. 



component of the electric field plays a role, by ensuring 
the validity of the condition (|T5|) . The z component of 
the magnetic fields does not affect the ray propagation, 
since it exerts a vanishing Lorentz force. However, the 
latter may affect the spin precession. If we assume that 
the input electrons have a spin parallel to the z axis, then 
the z component of the fringe magnetic field does noth- 
ing. However, the spin may be somewhat rotated away 
from the z axis by the transverse fringe magnetic field, 
and in this case the longitudinal component will induce 
some non-uniform precession around the z axis that may 
partly depolarize the beam. This effect will however be 
proportional to the time spent by the particles in the 
fringe field region, and hence will be negligible in the 
large length-to-gap ratio limit. 

We may also consider the idealized limit of an ex- 
tremely thin fringe field region in which the electric and 
magnetic fields vary rapidly from the vanishing values 
they have away from the Wien filter to the values they 
acquire well inside the filter, which then remain perfectly 
constant. This is the SCOFF limit mentioned above. In 
particular, in the SCOFF limit the vector potential A 
must go from zero to the value A — [0, 0, A z (x, y)] de- 
scribing the multipolar magnetic field inside the filter. 
This rapid transition does not give rise to any magnetic 
field, because it is curl-free. Therefore, the SCOFF limit 
has vanishing magnetic fringe fields and therefore abso- 
lutely no effect on the spin, nor on the lensing behavior 
of the filter. On the other hand, in the SCOFF limit the 
electric potential must go from zero to the nonzero value 
<fr(x,y) describing the multipolar electric field inside the 
filter. This in turn implies the presence of a very strong 
longitudinal electric fringe field that is responsible for the 
change of velocity given by Eq. (fT5|) . We stress that this 
idealized SCOFF description of the fringe fields needs not 
being realistic, but it is sufficient to identify the inelim- 
inable effects that the fringe field impose to the electrons, 
and therefore the possible fundamental limitations that 
may arise from them. 



Fringe field effects 

The actual fringe fields depend on the electrodes and 
magnetic poles geometry. It is however safe to assume 
that the fringe fields extend only for a longitudinal dis- 
tance along the z axis that is comparable to the trans- 
verse aperture of the Wien filter, and therefore much 
smaller than the filter length, in the large length-to-gap 
ratio limit. Moreover, the transverse components of the 
fringe fields can be taken to have exactly the same mul- 
tipolar geometry as the fields inside the Wien, and hence 
will give rise to similar effects both in lensing and spin 
precession. Their overall effect is therefore only to change 
the effective length of the filter. However, in addition 
there will also be a longitudinal z component that must 
be taken into account. In the lensing effect, only the z 



Layout of a spin-polarized TEM microscope 

We report in Fig. [5] the sketch of a possible TEM mi- 
croscope incorporating the proposed spin-polarization fil- 
ter based on the g-filter. Such a device could be used to 
perform spin-sensitive TEM microscopy on magnetic ma- 
terials or spintronic devices. It includes the four stages 
discussed in the main paper, i.e. (i) the fork hologram 
to set the initial OAM, (ii) the ^-filter to couple OAM 
to spin, (iii) an imaging stage to go in the Fourier plane, 
(iv) an aperture to select a single spinorbit component, 
as well as the additional electron optics needed for imag- 
ing and collimating the beam. The input aperture of the 
system is defined by the fork hologram mask and has a 
typical radius of few microns |2J . The output iris needed 
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FIG. 5. Sketch of a possible TEM microscope incorporating 
the proposed spin-polarization filter based on the g-filter. 

to select a single spin polarization can be set to a radius 
of few tens of microns, by carefully selecting the aperture 
plane position with respect to the focal plane. It will be 
presumably convenient to work slightly off focus to ob- 
tain larger beam waists without the need to introduce 
additional magnifying optics. 

We note here in passing that in the new generation 
of TEM microscopes with aberration-corrected probe a 
larger number of lenses are available for beam control, 
and in particular a pair of hexapoles. It could be worth 
then exploring the possibility of using directly these 
lenses as an effective g-filter, without introducing addi- 
tional optical elements. In this configuration, the mag- 
netic hexapoles with opposite sign could be used to pro- 
duce the polarization, while also compensating the intro- 
duced aberration. However stability issues might arise 
for large fields and the effect of fringing fields on the spin 
may also not be negligible, in this case, as the length- 
to-gap ratio would not be very large. Therefore, further 
analysis and simulations will be necessary in order to as- 
sess this possibility. 



The role of electron diffraction 

In our treatment we fully included the diffraction ef- 
fects only in the free propagation (through a lensing sys- 
tem) taking place after the g-filter device. Indeed, this 
stage is where the correlation between OAM and radial 
profile develops and is then exploited for separating the 
different spin components. This free propagation, end- 
ing in the far-field and giving rise to the beam profiles 



shown in Fig. 3 of the main paper, has been evaluated 
using a fully quantum mechanical treatment. The diffrac- 
tion effects are obviously important in the propagation 
from near to far field, and in particular in determining 
the presence or absence of destructive interference effects 
at the vortex position that are essential in our proposed 
scheme. 

One may wonder if we can legitimately ignore the 
diffraction effects in the other stages of our proposed 
setup. For the propagation through the q-filter, we could 
assess directly the role of diffraction in the case of homo- 
geneous field, for which we have both the exact solution 
of Pauli's equation (including diffraction) and the classi- 
cal trajectory solutions. We find that diffraction is only 
relevant if we look at the electron beam very close to a 
focal point, where its width becomes of the order of few 
nanometers. Away from the focus, the ensemble of clas- 
sical trajectories reproduces perfectly well the quantum 
beam evolution. This result is not surprising, as in the q- 
filter all length scales are 8-9 orders of magnitude larger 
than the electron wavelength A = 27r/fc c , which is of few 
picometers for energies of the order of 100 keV (typical of 
electron microscopes). Therefore, we may safely assume 
that diffraction can be ignored also in inhomogeneous q- 
filter geometries, as long as no focusing occurs inside the 
filter. Our classical ray-tracing simulations confirm that 
we avoid a focusing inside the g-filter. It should be also 
noted that in designing and analyzing electron-optics el- 
ements as those commonly used in electron microscopy, 
it is quite standard to ignore quantum effects and use 
classical ray tracing or geometrical optics. 

The input and output apertures, being again about 7 
orders of magnitude larger than the electron wavelength, 
will also give negligible diffraction effects (only some ra- 
dial fringes extremely close to the aperture edges are ex- 
pected, which should not affect the proper working of the 
proposed setup). 

Violation of Bohr's "theorem" 

As reported in Ref. [13], few years after the discov- 
ery of electron spin, several experiments were proposed 
(and sometimes attempted) aimed at separating electrons 
according to their spin, mainly by exploiting inhomoge- 
neous magnetic fields acting on the magnetic dipole asso- 
ciated with it (similar to the Stern-Gerlach experiment). 
In 1929, Niels Bohr demonstrated a sort of "impossibility 
theorem" [13], stating that no Stern- Gerlach-like exper- 
iment (including Knauer's proposal, in which an elec- 
tric field was used to balance the Lorentz force) could be 
used to separate electrons according to their spin by more 
than a small fraction of the beam transverse extension. 
This was related to the uncertainty principle and to the 
unavoidable orthogonal magnetic field accompanying the 
required magnetic field gradient. A simple explanation 
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of this result is that the spin magnetic moment vanishes 
in the classical limit fi — ¥ 0, so all attempts at separating 
particles according to the magnetic-moment force occur- 
ring in a magnetic field gradient give rise to a vanishing 
separation in the classical limit, and hence a separation 
that is comparable to quantum uncertainty. Actually, 
this statement has been later proved not to be strictly 
valid (see, e.g., Refs. [H El), 

although it remains ex- 
tremely hard to overcome in practice. 

Our approach is however fundamentally different from 
past attempts. Indeed, we exploit the quantum nature 
itself of the electrons to operate the spin separation (as 
such, our proposal bears some similarities with a very 
recently proposed scheme for spin-polarization in which 
electron diffraction is also exploited 20]). As we have 
seen, the magnetic dipole force in our device acts only in 
the azimuthal direction, which is not the direction along 
which we separate the electrons. The torque associated 
with this force imparts the OAM variation correspond- 
ing to the azimuthal geometric phase. However, at the 
exit of our Wien filter the electrons show no significant 
spatial separation according to the spin. It is then the 
subsequent free propagation and diffraction that acts to 
separate the electrons along the radial direction, which 
is not the direction of the main field gradient. In par- 
ticular, an efficient radial separation is ensured by the 
presence of a vortex at the beam axis for one spin com- 
ponent and not for the other. In other words, we use the 
quantum wavy nature of the electrons to separate them 
according to the spin, and hence our method is not sub- 
ject to Bohr's impossibility theorem. As a further argu- 
ment, even in the classical Stern-Gerlach approach Bohr's 
theorem does not forbid the creation of spin-polarized 
components at the edges of the beam, so that by suit- 
able aperturing it should be possible to generate spin- 
polarized beams, although with very low efficiency. Our 
method is also based on selecting only an edge of the 
beam, although this is done in the radial direction and 
the efficiency is strongly improved by the vortex-related 



destructive interference effects. Finally, we note that in 
our geometry the unavoidable orthogonal magnetic field 
associated with the azimuthal field gradient generating 
the dipole forces, which is the main source of problems 
in the Stern-Gerlach approach to spin-polarized electron 
separation, is already taken into account and corresponds 
to the radial magnetic field of the quadrupole or hexapole 
configuration. So its effect has been already fully consid- 
ered in our treatment. 
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